Robust quantum control by shaped pulse 
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Considering the problem of the control of a two-state quantum system by an external field, we 
establish a general and versatile method that allows the derivation of smooth pulses, suitable for 
ultrafast applications, that feature the properties of high-fidelity, robustness, and low area. Such 
shaped pulses can be viewed as a single-shot generalization of the composite pulse-sequence technique 
with a time-dependent phase. 

PACS numbers: 03.65.Aa, 32.80.Qk, 42.50.Dv 42.50.Ex 



o 

(N 

< 



Or 
i , 

H— > , 

a: 

cr 



> 

o 
o 



* 



Modern applications of quantum control, such as quan- 
tum information processing [l|, require time-dependent 
schemes featuring three important issues: The transfer 
to the target state should be achieved (i) with a high fi- 
delity, typically with an admissible error lower than 10~ 4 
l], (ii) in a robust way with respect to the imperfect 
knowledge of the system or to variations in experimental 
parameters, and (iii) with a minimum time of interac- 
tion and low field energy in order to prevent unwanted 
destructive intensity effects. 

The Rabi method (see e.g. Q), corresponding to an ex- 
act resonant coupling between two quantum states, leads 
to a complete transfer to an excited state from a ground 
state by a Rabi frequency of area A = n (7r-pulse tech- 
nique). This defines the transfer quantum speed limit in 
the sense that the 7r-area of the Rabi frequency is the 
smallest area that gives a complete transfer P = 1 Q. 
Any (time-dependent) variation in the detuning requires 
a Rabi frequency area larger than 7r, i.e. an increase of 
the fluence of the field, in order to recover the complete 
transfer. This extra energy can be used to satisfy ad- 
ditional constraints such as robustness with respect to 
the variations of parameters. In this framework, adia- 
batic techniques are famous examples [J]; they however 
require in principle a large pulse area and do not lead 
to an exact transfer. Improvements to optimize its effi- 
ciency by parallel adiabatic passage [5| or by shortcuts 
to adiabaticity [6] have been proposed. Optimized shap- 
ing of the detuning and of the field amplitude leading 
to both robust and precise transfer remains however an 
open question. 

Making the transfer robust means that the errors that 
can appear in control parameters have to be compen- 
sated, without even knowing these errors. A practi- 
cal measure of the robustness can be defined through 
the deviation of the excitation profile as a function of 
the considered parameters. The use of composite pulses 
[j| is a popular method for self-compensation of errors. 
The basic idea consists in replacing the single resonant 
pulse driving the quantum transition by a sequence of 
pulses with well-defined static phases. These phases, as 



* sguerin@u-bourgognc.fr 



well as the area of the individual pulses, can be deter- 
mined so that the derivatives of the excitation profile 
are nullified order by order by increase of the number 
of pulses. Recently, an SU(2) algebraic approach has 
been extended for the design of composite sequences of 
pulses with smooth temporal shapes and constant 8] or 
time-dependent Q detuning. We emphasize that such 
smooth pulses are more appropriate for ultrafast pro- 
cesses, for instance, with femtosecond pulses, requiring 
their production in the spectral domain [1J| , since rectan- 
gular pulses would lead to a prohibitively large spectrum. 
Typically, for a resonant complete inversion within a two- 
state system, while the profile P deviates around A = it 
as P ~ 1 — (A — it) 2 with respect to the area for the Rabi 
method, it deviates around A = 7r as P ~ 1 — (A — ix) 2n 
with n the (odd) number of composite resonant n— pulses. 
Techniques of optimal control are also actively developed 
for this purpose of robustness [11| , but they only lead to 
some purely numerical solutions without any insight into 
the physical mechanism nor the analytical derivation of 
the control fields. 

In this Letter, we establish a control strategy that al- 
lows a robust transfer by a pulse specifically shaped in 
phase and amplitude. This can be viewed as a general- 
ization of the strategy relying on a series of composite 
pulses with static phases to a single-shot pulse of time- 
dependent phase. We first show that the issue of robust- 
ness can be reduced to nullifying the derivatives of the 
excitation profile to a target state order by order. The 
central result of this work is that this can be in general 
achieved by an oscillatory parameterization of the phase 
of the wave function. This continuous trigonometric ba- 
sis is the key difference with respect to known methods 
which make use of stepwise functions (e.g. for the phase 
in the case of composite pulses) . This allows an explicit 
derivation of the components of the shaped pulse. Our 
technique is explicitly shown for the robust complete in- 
version with respect to the pulse area, to the detuning, or 
to both parameters. This approach is however versatile 
and can be be applied to other types of robustness and 
to the transfer to more complicated targets, such as fully 
robust quantum gates in view of applications in quantum 
information processing. 

The resulting smooth shaped pulses feature the re- 
quired properties of high-fidelity, high-order robustness, 



and low area, in particular smaller than their compos- 
ite counterpart for the same efficiency. Furthermore the 
pulses we derive have an explicit analytic form with very 
few parameters to adjust contrary to numerical optimal 
control procedures where a large number of parameters 
is used (see for instance [12|]), and the resulting solution 
features a much simpler form than the ones usually ob- 
tained diii!]. 

Our discussion is based on a resonant system (rotat- 
ing wave approximation) between two states |1) and |2), 
for which the most general Hamiltonian governing the 
dynamics can be written [l5[ 



TDSE, we get: 
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The Rabi frequency is decomposed into an absolute value 
Q,(t) > (proportional to the field amplitude for a one- 
photon transition), of area A = J Q(t)dt, and a phase 
T](t). The detuning A(i) between the field frequency and 
the transition features a time- independent (static) de- 
tuning and a time-dependent part, produced, at the fem- 
tosecond timescale, by a spectral phase shaping through 
the phase r\ [10l | . 

The errors arise from the imperfect knowledge of the 
area A (for instance through an imperfect knowledge of 
the coupling constant), of the static detuning 6 (often re- 
ferring to a problem with an inhomogeneous broadening 
of an ensemble due to its environment), or to dynam- 
ical fluctuations of the pulse shape or of its instanta- 
neous phase (corresponding to fluctuations of the time- 
dependent part of the detuning). For simplicity, we focus 
our discussion on the deviation with respect to the area A 
of the Rabi frequency and the static detuning S. Robust 
methods are designed to improve the quadratic deviation 
benchmark of the Rabi method. 

The solution of the time dependent Schrodinger equa- 
tion (TDSE) ih£r<f> = H(j> can be parameterized in 
the most general way with two angles {8 = 8(t) G 
[0, 7r], ip — 77 = <p(t) — i](t) e] — 7r, 7r]} on the Bloch sphere 
of cartesian coordinates (with p mn = (m\ip) (ip\n)) p x — 
P21 + P12 = sin 9 cos ip, p y = i(p 2 i — P12) = sin 9 sin ip, 
Pz = P11 — P22 = cos#, and with a global phase 7 = 7(t) 
as 
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We first 
diag[e^"'/ 2 , 
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a phase transformation T = 
to deal with a real symmetric 



T^HT - ihT^ 



dT 
dt 



H(A + rj,n,0) 



Hamiltonian H = 

of solution ip = T^4>. This shows that the phase ij 
can be incorporated in the detuning term and simply 
interpreted as the rotation of the axes x and y about the 
z-axis of the Bloch sphere. One can thus assume r\ = 
without loss of generality. Inserting this solution in the 



9 = fl sin (/?, 

ip = A + ft cos ip cotan 9, 
cos ip ■ cotan ip 
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(3a) 
(3b) 

(3c) 



We assume <fi(ti) = 1) as initial condition, correspond- 
ing to the initial conditions 0(U) = 9 h = (north 
pole), p(U) = p, = j(ti) = 74 (not specified by 
the initial state). We consider the complete popula- 
tion transfer, which is achieved for the final condition 
9(tf) = Of = it (south pole). For any trajectory fea- 
tured by <p(t), <r- 0(t) — » 7r on the Bloch sphere, one 



can integrate ([3a] 



9f — 0i = 7r = J > ds Ct(s) sin tp(s). 

Since J f dsil(s) sin ip(s) < f.f ds£l(s), one concludes 
that a complete population transfer is achieved when 
J. f dsfl(s) > it, in consistency with Ref. [3j. The min- 
imum area it is obtained for the meridian p = 7r/2, im- 
plying 7j = 7r/2, 7 = from (f3"c|) and A = from (|3b| . 
which corresponds to the Rabi transfer. On the other 
hand, the ideal adiabatic solution is derived when A ^> 1 
giving ip — > 0,7r [in order to have a bounded from pap ], 
that is 7(t) ->■ ±fl ^/fi(s) 2 + A(s) 2 ds from ([5c]). From 
@, we recover the well-known adiabatic dynamics. 

More generally, the transfer to a given target state up 
to a (global) phase corresponds to the desired final con- 
ditions 9 f and ipf . The control of the phase of the target 
state would require the additional condition j(tf) = 7/, 
which will not be considered here. We will additionally 
assume Cl(ti) = fl(tf) — in view of a practical imple- 
mentation. We will also consider a monotonic increasing 
of 9, i.e. 8 > 0, with ft > 0, which leads to < ip < it 
from Eq. ([3all. 

The goal of the control consists in finding a particu- 
lar solution of © with the final conditions correspond- 
ing to the target state, which is robust and of lowest f2 
area. The strategy is similar to the one used for generat- 
ing composite pulses: One nullifies the derivatives of the 
transfer profile with respect to the considered parame- 
ters. This can be here achieved with the Hamiltonian of 
the form 



H*At) = 2 



-A(i) Q(t) 

n(t) A(t) 



-S afl(t) 
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The corresponding solution 4>a,s(t) is parameterized by 
a and 5 for given functions fl(t) and A(t). We assume 
that fl(t) and A(t) are such that the complete transfer 
to the target state 4>t is achieved at t = tf for S = 
and a — 0: <f>o(tf) = 4t- O ne makes a perturbative 
expansion of 4> a ,s(tf) with respect to a and S taking the 
second matrix term of (UJ) as a perturbation denoted V 
(the small parameters are both a and S assumed of the 
same order) \v^: (<h\<t>a,s(tf)) = 1 + °i + 2 + 3 + - ■ ■ , 
where O n denotes the term of total order n, giving for 
the excitation profile 



\(<h\<l> a ,s(fif))\ 2 = 1 + Ox + o 2 + o 3 + 



(5) 



with O n the term of order n. The first two terms read 

O x = -i [ f (MtW(t)\Mt))dt = -i [ ! e(t)dt, (6a) 
Jti Ju 

2 = (-i) 2 f'dtf dt'[e(t)e(t') + f(t)f(t% (6b) 



r , (7) 



ti ju 



with 0\ = 0± + 0\, e = —h (8 cos 9 — 07 sin 2 



/=(^o|F|^) = - 



<5 sin ( 



ai —7 sin 20 — 19 



and the orthogonal solution of the TDSE 
(j)x(t) = [e^/ 2 sin(6>/2),- e -^/ 2 cos(6>/2)] T e^/ 2 such 
that (4>±(t)\<j)o(t)) = 0. The other terms can be deter- 
mined from the symbolic diagrams depicted in Fig. [T] 
Since e(t) is real, there is no first-order deviation for the 
excitation profile. One can simplify the second order as 



2 = 02 + 2 + O x Ox = - 



f(t)dt 



(8) 



£dtj*dt'[a(t)b(t') + 

This property 



property 
£a(t)dtf?b(t)dt 
J^dt a(t)f^dtb(t)f*dt'c{t' 



using the 
a(t')b(t)] = 
also implies 

/ T T dt f* df f r '' dt" [a(t)b(t')c(t") + 

b(t)a(t')c(t") T + b(t)c(t')a(t")} and extends as 

jf dt f* dt' £ dt" J2 a a(t( ff (°)))6(i( CT ( 1 ))) C (i( CT ( 2 ))) 

/ a(t)dtf b(t)dt J c{t)dt, where J2 a means the 
summation over the six permutations of the number of 
primes (between and 2). This is used to determine 
relatively simple integrals for higher orders. The third 
order reads: 



if 



dt / dt' / dt"Im[f(t)e(t')f(t")}. (9) 



O, 



Robustness at a given order n is obtained when the pa- 
rameters of the field are chosen such that they allow nul- 
lifying the integrals O m for m < n. The second-order 
robustness issue corresponds to the two equations: 



e^ sm9 dt = 0, 



7 „ 



ife l7 sin26ll 



for the robustness with respect to the detuning 5, and to 
the pulse area respectively. If one considers the robust- 
ness only with respect to the pulse area, the correspond- 
ing equation involves an integral which is independent of 
the particular temporal parameterization of 9. For some 
integrals, such as, for instance the second order robust- 
ness with respect to the detuning [left equation of (TTU)) ]. 
we need additionally an explicit time parameterization of 
9. For the Rabi method (A = 0), these equations (flQ)) 
cannot be satisfied for any pulse shape: It is robust up 
to the first order for the excitation profile. Below, we ex- 
plicitely derive solutions for the issue of inversion (with 
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FIG. 1. Symbolic path-diagrams giving the construction of 
the On integrals. The symbol e stands for e(t), e' for e(t'), 
and so on. For instance, for n — 3 (n = 4), the diagram 
features 4 (8) paths. Its extension for larger n is direct. 



9 varying from to 7r). The same technique applies for 
other target states. 

These equations (flU)) (and the ones corresponding to 
higher orders) can be exactly solved by an oscillatory 
parameterization with a Fourier series of the global phase 
as a function of 9, 7(4) = j(9) (with the label a or b to 
identify it): 



7 a (6») = 29 + d sin(20) 
76 (6») = 9 + d sin(20) 4 



f C„ sin(2n0) 4 
C„ sin(2n0) + 



(Ha) 
(lib) 



choosing the smooth parameterization 9{t) = 
7r(erf(i/T) + l)/2 (which will lead to a smooth 
pulse) with T featuring the duration of interaction. 
These parameterizations (|11[) allow a simple expression 
in view of controlling the phases, since, from ([5c]). we 
have cotant^ = sin(9)cFf/d6. These choices give 7/ = 2ir 
(it) for parameterization (jllap [(jTTbJ)]. For the choice 
(lllal) . the imaginary parts of Eq. (JTUJ) are nullified for 
all C„'s. On the other hand, for (|llbl) . the real parts 
are nullified for all C n 's. In both cases, one has then to 
adjust the coefficients C n to nullify the other integrals. 
An important result is that, due to the oscillatory nature 
of the parameterization, only a few coefficients are 
needed to nullify the integrals; typically one additional 
coefficient different from zero is considered to nullify 
another single (real) integral. 



TABLE I. Robustness of order n (i.e. O m < n = 0) with respect 
to the area (referred to as type A), to the detuning (type S), 
or both (type AS), with the parameterization (jllap (referred 
to as param. a) or (|llb[l (param. b), and the coefficients Cj, 
j — 1, 2, 3, (Cj>3 = 0). We have considered parameterizations 
leading to low pulse areas. 



Type Param. 


Order 


d 


c 2 


C3 Pulse area (X7r) 


A a 


3 


-1 





2.16 


A b 


3 


-1.6788 





2.09 


5 a 


3 


-0.2305 





1.78 


AS b 


2 


-1.189 


0.7285 


2.23 


A a 


5 


-2.4864 


-0.74 


3.14 


A a 


7 


-3.46 


-1.365 


-0.5 3.86 



Some obtained coefficients are summarized in 
Table Q] The resulting pulse area J t f £l(t)dt = 



JT yl + (pf) 2 sin 2 dO is mentioned. One can no- 
tice increasing areas when more coefficients different 
from zero are considered. Table U shows for instance 
that one obtains a pulse area of only 1.787T for the third 
order robustness [corresponding to the integrals © 
and © nullified, i.e. up to the third order] solely with 
respect to the detuning taking the parameterization 



(jllap with C\ = —0.2305. Alternatively, we obtain at 
best a pulse area of 2.09-7T for the third order robustness 
solely with respect to the pulse area taking the param- 
eterization (lllbp . We emphasize that this latter choice 
of parameterization allows a smooth pulse of slightly 
smaller area than the non-smooth pulse proposed in [16j 
(with C\ = —1, see first line of Table [J). In both cases, 
the third order integrals ([9]) are systematically nullified 
when the coefficients are adjusted to nullify the second 
order ((5]) due to the symmetry of the parameterization. 
One can alternatively force robustness with respect to 
both the detuning 5 and the pulse area, nullifying both 
terms of Eq. (JTUJ) (see fourth line of Table fl} . This shows 
the remarkable result that the obtained pulse area is 
only slightly larger than the one obtained above for 
robustness solely with respect to the pulse area (but for 
a robustness of order 2). 

If one considers robustness at high order, one has to 
nullify the higher order terms. The results for the ro- 
bustness with respect to the pulse area at orders 5 and 7 
are shown in Table HI 
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FIG. 2. Robustness with respect to the pulse area: Popu- 
lation transfer P2 at the end of the pulse as a function of 
the relative deviation a (dimensionless) of the area. Different 
orders of robustness (from order 3 to 7, as indicated next to 
the curves), where the used coefficients and parameterizations 
are indicated in Table [H are compared to the standard Rabi 
method. Upper: Linear scale of P2 ■ Lower: Logarithm to the 
basis 10 of the deviation. 



The robustness with respect to pulse area at different 
orders is demonstrated in Fig. [51 It shows the numerical 
deviation of the transfer profile for various values of a 
as defined in Eq. (U) (here we take S — 0), with the pa- 
rameters of Table U and the parameterization (jllap . As 
expected, the transfer profile becomes flatter and flatter 
when one nullifies higher order terms in the perturba- 
tive expansion (|S|) (see upper frame of Fig. [2]). We can 
notice that the profile is not symmetric with respect to 
a = 0. More precisely, the transfer is less robust for 
a smaller area (corresponding to negative values of a). 
This result is in accordance with the expected property 
that the robustness is better for a larger pulse area. The 
lower frame Fig. [5] depicts the deviation of the excitation 
profile at a logarithmic scale. It shows the remarkable 
result that the population transfer is accomplished with 
the 10~ 4 high-fidelity accuracy benchmark (indicated by 
the horizontal dashed line in Fig. 0) even with an error 
in area up to 17% for the highest order solution (limited 
by a negative area deviation). 




-0.5 0.5 
t (in unit of T) 

FIG. 3. Derived Rabi frequency pulse shapes for robustness 
of orders 3, 5 and 7 with respect to the pulse area (with the 
parameters of Table [I]) . 



Figure [3] displays the shapes of the Rabi frequencies 
that are derived: They oscillate more and more with a 
larger area for a higher order robustness. These oscil- 
lating pulses can be interpreted as an efficient general- 
ization of the composite pulse sequence, but with time- 
dependent phases instead of the static ones used for the 
composite pulses. 

In conclusion, we have constructed and analyzed a 
generalization of the technique of composite pulse se- 
quence with static phases to a single-shot pulse of time- 
dependent phase. We have reduced the robustness is- 
sue to a problem of nullifying integrals, and have shown 
that this can be achieved for an oscillatory parameteri- 
zation of the phase of the solution. The resulting smooth 
shaped pulse, suitable for ultrafast applications, feature 
properties of high-fidelity, high-order robustness, and low 
area. Furthermore the pulses we derive have an explicit 



and relatively simple form which could be easily imple- 
mented experimentally. A complete analysis has been 
given for the robustness with respect to the pulse area, 
to the detuning, or to both parameters, and an explicit 
parameterization has been demonstrated for the inver- 
sion problem. Here only two or three free parameters C n 
have been systematically adjusted to achieve a very effi- 
cient robust inversion. The techniques which we present 
is however versatile and can be applied to other types of 
robustness and target, such as coherent superposition of 
states or simultaneous control. Using the same param- 
eterization, one can for instance derive a shaped pulse, 
with a few coefficients to be adjusted, producing a robust 
superposition of state [with 6 varying from to a given 
value featuring the weight of the desired superposition, 
and iff chosen from (|3c[) ] . Ultimately, using the same 
technique (but adapting the parameterization), one can 
also produce robust propagators with self-compensation 



of errors, suitable for quantum information processing. 
In this case, more integrals to be nullified will be in- 
volved (but of the same type that the ones described in 
this work) , which will need more coefficients C„ to be ad- 
justed: algorithms of optimal control will be then needed 
to optimize these coefficients [17]. This task will however 
be relatively easy in comparison with traditional optimal 
control techniques which need in general more than hun- 
dred coefficients to be adjusted. This simplification will 
be also useful in view of extension to more complicated 
systems. 
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